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Within the framework of unified approach we study the Casimir-Lifshitz interaction, the van der 
Waals friction force and the radiative heat transfer at nonequihbrium conditions, when the inter- 
acting bodies are at different temperatures, and they move relative to each other with the arbitrary 
velocity V. The analysis is focused on the surface-surface and surface-particle configuration. Our 
approach is based on the exact solution of electromagnetic problem about the determination of 
^^ ' the fiuctuating electromagnetic field in the vacuum gap between two flat parallel surfaces moving 

Cn I relative to each other with the arbitrary velocity V. The velocity dependence of the considered phe- 

^.^. nomena is determined by Doppler shift and can be strong for resonant photon tunneling between 

j^ ' surface modes. We show that relativistic effects give rise to a mixing of the contributions from the 

electromagnetic waves with different polarization to the heat transfer and the interaction forces. We 
find that these effects are of the order (V/c)^. The limiting case when one of the bodies is sufficiently 
w-v ' rarefied gives the heat transfer and the interaction forces between a moving small particle and a 

surface. We also calculate the friction force acting on a particle moving with an arbitrary velocity 
relative to the black body radiation. 
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I. INTRODUCTION 



d • All bodies are surrounded by a fluctuating electromagnetic field due to the thermal and quantum fluctuations of the 

charge and current density inside the bodies. Outside the bodies this fluctuating electromagnetic field exists partly 

T ■ in the form of propagating electromagnetic waves and partly in the form of evanescent waves. The theory of the 

' O ' fluctuating electromagnetic field was developed by Ryto\i"— . A great variety of phenomena such as Casimir-Lifshitz 

5 , forces^, near-field radiative heat transfer—, and friction forces^ can be described using this theory. 

Lifshita^ used the Rytov's theory to formulate a very general theory of the dispersion interaction in the framework 
of the statistical physics and macroscopic electrodynamics. The Lifshitz theory provides a common tool to deal with 
dispersive forces in different field of science (physics, biology, chemistry) and technology. 
^^ ' The Lifshitz theory is formulated for systems at thermal equilibrium. At present there is an interest in the study 

I of systems out of the thermal equilibrium (see^ and reference therein) . The principal interest in the study of systems 

r-^ out of thermal equilibrium is connected to the possibility of tuning the interaction in both strength and sign*'^. Such 
("^ , systems also give a way to explore the role of thermal fiuctuations usually masked at thermal equilibrium by the 
T-H ' T = K component which dominates the interaction up to very large distances, where the actual total force result in 

f I , to be very small. The Casimir-Lifshitz force was measured at very large distances and it was shown that the thermal 

f~^ ' effects of the Casimir-Lifshitz interaction are in agreement with the theoretical prediction^. This measurement was 
QQ ] done out of thermal equilibrium, where thermal effects are stronger. 

^^ . Further thermal non-equilibrium effects were explored by Polder and Van Hove^, who calculated the heat-flux 

l~^ ' between two parallel plates. At present there is an increasing interest for studying near-fleld radiative heat transfer—"— 

• i-H , which is connected with the development of a near-field scanning thermal microscope^. The existing studies are 

/\^ • limited mostly by the case when the interacting bodies have different temperatures but they are at rest. For recent 

J^ , review of near-fleld radiative heat transfer between bodies, which are at rest, se o"'"^'^^ . 

■ - - ' Another non-equilibrium effects are realized for bodies moving relative to each other. In^ we used the dynamical 

modification of the Lifshitz theory to calculate the friction force between two plane parallel surfaces in relative motion 

with velocity V. The calculation of the van der Waals friction is more complicated than of the Casimir-Lifshitz force 

and the radiative heat transfer because it requires the determination of the electromagnetic field between moving 

boundaries. The solution can be found by writing the boundary conditions on the surface of each body in the rest 

reference frame of this body. The relation between the electromagnetic fields in the different reference frames is 

determined by the Lorenz transformation. In^ the electromagnetic field in the vacuum gap between the bodies was 

calculated to linear order in V/c. It was shown that linear terms in the electromagnetic field give the contribution to 

the friction force of the order (V/c)'^. Thus, these linear terms were neglected inS and the resulting formula for friction 
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FIG. 1: Two semi-infinite bodies with plane parallel surfaces separated by a distance d. 
other with velocity V . 



The upper solids moves parallel to 



force is accurate to order (V/c)'^. The same approximation was used inlS to calculate the frictional drag between 
quantum wells, and in2Si2i to calculate the friction force between plane parallel surfaces in normal relative motion. 
For a recent review of the van der Waals friction see^^. 

In this paper within the framework of unified approach we study the Casimir-Lifshitz interaction, the van der 
Waals friction force and the radiative heat transfer at nonequilibrium conditions, when the interacting bodies are 
at different temperatures, and they move relative to each other with the arbitrary velocity V. Our study is focused 
on the surface-surface and surface-particle configuration. In comparison with previous studies we consider more 
general nonequilibrium conditions. In the existing literature the Casimir-Lifshitz interaction and the radiative heat 
transfer for the surface-surface configuration were studied only for the systems out of the thermal equilibrium^ " ^^' ^^ . 
The van der Waals friction is studied for this configuration only for systems at thermal equilibrium^^. In Sec. Hi] 
we calculate the fluctuating electromagnetic field in the vacuum gap between two plane parallel surfaces, moving 
in parallel relative to each other with arbitrary velocity V. In comparison with previous calculations^ii2ri2i we do 
not make any approximation in the Lorentz transformation of the electromagnetic field by means of which we can 
determine the field in one inertial reference frame, knowing the same field in another reference frame. Thus, our 
solution of the electromagnetic problem is exact. Knowing the electromagnetic field we calculate the stress tensor 
and the Poynting vector which determine the interaction force and the heat transfer, respectively. We calculate the 
friction force and the conservative force, and the radiative heat transfer in Sects. IIIIl IIVI and IVl respectively. Upon 
going to the limit when one of the bodies is rarefied we obtain the interaction force and the heat transfer for a small 
particle-surface configuration. In Sec. IVII we calculate the friction force on a small particle moving relative to black 
body radiation. The same problem was considered in^^. In comparison with this study our treatment is relativistic 
and we take into account the contribution not only from the electric dipole moment but also from the magnetic 
moment of the particle. Recently^- it was shown that the magnetic moment gives the most important contribution 
to the near-field radiative heat transfer for metallic particles. The same is true for the friction force. The conclusions 
are given in Sec. IVIII 



II. CALCULATION OF THE FLUCTUATING ELECTROMAGNETIC FIELD 



We consider two semi-infinite solids having flat parallel surfaces separated by a distance d and moving with velocity 
V relative to each other, see Fig. [TJ We introduce the two coordinate systems K and K' with coordinate axes xyz 
and x'y'z' . In the K system body 1 is at rest while body 2 is moving with the velocity V along the x— axis {xy and 
x'y' planes are in the surface of body 1, x and x'- axes have the same direction, and the z and z'— axes point toward 
body 2). In the K' system body 2 is at rest while body 1 is moving with velocity —V along the x— axis. Since the 
system is translational invariant in the x = (a;, y) plane, the electromagnetic fleld can be represented by the Fourier 



integrals 

E(x, z, t) = J^ dw J ^e^'i— *E(q, cu, z), (1) 

B{^,z,t)^ p d^.y^e^'i— *B(q,c.,z), (2) 

where E and B are the electric and magnetic induction field, respectively, and q is the two-dimensional wave vector 
in xy- plane. After Fourier transformation it is convenient to decompose the electromagnetic field into s- and p- 
polarized components. For the p- and s -polarized electromagnetic waves the electric field E(q, a;,2;) is in plane of 
incidence, and perpendicular to that plane, respectively. In the vacuum gap between the bodies the electric field 
E(q, a;,z), and the magnetic induction field B(q, w,2) can be written in the form 

E(q, uj, z) = (vshs + Vpn+) e'''^' + {wsfis + Wpfip) e"*''^^ (3) 



B(q,w,z) = - [([k+ X n,]vs + [k+ x n+]vp) e'' 



+ ([k X ns]ws + [k xfip ]wp) e ''''^] (4) 

where k± = q±zk^, k^ = {{uj/cf-q'^y/^, n^ = [zxq] = (-9y,fe,0)/q, fi^ = [fc^^xn^] = {Tqxkz,Tqykz,q^)/{kq), k = 
uj/c. At the surfaces of the bodies the amplitude of the outgoing electromagnetic wave must be equal to the amplitude 
of the reflected wave plus the amplitude of the radiated wave. Thus, the boundary conditions for the electromagnetic 
field at z = in the K- reference frame can be written in the form 

Vp(s) = Ripis) (w, q)wp^s) + ^fp(^) {^, q) (5) 

where RipU) i'-^) is the refiection amplitude for surface 1 for the p{s) - polarized electromagnetic field, and where 
El u) {^) is the amplitude of the fluctuating electric field radiated by body 1 for a p(s)-polarized wave. In the K - 
reference frame the electric field can be written in the form 

E'(q',c.', z) = {v',n', + v'^h'+) e^^^' + {w',n', + w'^n'^) e"^^'^ (6) 



where q' = {q'^,qy,Q), g^ = {q^ - ^fc)7, uj' ^ {uj - Ffe)7, 7 = 1/^1-/32, j3 = V/c, h'^ = {-qy,q'^,0)/q', n'^ 
{Tq'xkz,Tqykz,q''^)/{k'q'), 



q' - 7y^92-2/3fcfe + /32(P-g2)_ 
The boundary conditions at z = d in the K - reference frame can be written in a form similar to Eq. ^: 



2ik^d-Q^ , l,,,l J\„J iJk^dT^lf /,/ / 



.p(,) - e--"i?2p(.)(^',g')<(.) + e^"^'*i?;;(,)(a.',g'), (7) 

where i?2p(s) {^) is the refiection amplitude for surface 2 for p{s) - polarized electromagnetic field, and where E^ ,^. {ui) 
is the amplitude of the fluctuating electric field radiated by body 2 for ap(s)-polarized wave. A Lorentz transformation 
for the electric field gives 

E'^ = E,, 4 = {Ey - (3Bzh, E', = {Ez + l3Byh (8) 

Using Eqs. pilSl) and © we get 

, fc'7 



V. 



J [-(3kzqyVs + {q^ - I3kqx)vp] , (9) 



^ kqq 

"'p = Y^ [P^-'iy^^ + (9' - Pkq,,)wp\ , (10) 

v's = T— 7 [PkzqyVp + {q^ - !ikq.^)Vs\ , (11) 



k''y 
w's = T — 7 [-Pk^qyWp + {q^ - l3kqx)ws 

Substituting Eqs. PIT2|) in Eq. and using Eq. (O we get 

(g^ - l3kqx)DppWp + /3k^qyDspWs 



(12) 



2ik^d r>l 



(g^ - fikq,)Elp - pk^qyEi 






(13) 



where 



(g^ - /3kqx)DssWs - jSk^qyOpaWp 



^2rk.dj^, 



{q^^f3kq^)El+l3k,qyEl 



fc'7 



2s ! 



-Dpp — 1 — e * ^ i?ipi?2p, £'ss = 1 — e * ^ RisR2s: 



Dsp — 1 + e ' ^ i?isi?2„, i'ps — 1 + e ' ^ i?ipi?2( 



-^2p(s) = -R2p(s)(w',g')- From Eqs. (|13I14I) and © we get 



(14) 



w;p = { [iq^ - pkq^fR'^pD,, - p^kXRLDsp] E^e 



f „2ik^d 



I \ Tpf ^2ik^d 



-PKqy{q^ ~ l3kq^){R'2p + RDELe 



kqq' r 

k'j 



{q^ - /3kq^)D,sE!/p - iSk^qyDspE'i] e^'^^^JA"! 



(15) 



Vp = { [{q^ - pkq^fD.s + P^k^yD^p] El 



-f3k,qy{q^ - pkq,)R,p{R',p + R'2,)e^'^'^Ei 



kqq' 

k'-f 



R 



ip 



(g2 _ l3kq,)D,sEg - l3k,qyDspEg 



Ak^d 



}A- 



(16) 



{ [(g2 _ pkq^fR'^.Dpp - P^klq^R'^pOp^] El^e^^"^'' 



+Pk,qy{q^ ~ Pkq,){R'2p + R'2,)Ele 



I \ 771/ „2ik^d 



kqq' 

k'"f 



(g2 - pkq,)DppE'l + pk^qyOp^E'^^^ e^'^^'^JA"!, 



(17) 



Vs = { [{q^ - Pkq^fDpp + fi^klqlDp,] E(^ 



+Pk,qy{q^ - pkq^)R,p{R'2p + R'2^)e^'^'''El 



kqq' 



where 



^i?i, [{q^ - f3kq^)DppEg + iSk^qyDpsEg] e^'^'-'^JA"! 



A^{q^- /3kq,)^D,,Dpp + P^klqlDp,D,p. 



(18) 



The fundamental characteristic of the fluctuating electromagnetic field is the correlation function, determining the 
average product of amplitudes £'■', >(q, oj). According to the general theory of the fluctuating electromagnetic field 
(see for a example^): 



+ (fc,-A::)(i?;(,)-i?p(,))] 



(19) 



where < ... > denote statistical average over the random field. We note that k^ is real for q < uj/c (propagating 
waves), and purely imaginary for q > uj/c (evanescent waves). The the Bose-Einstein factor 



n(w) 



1 



^hu/ksT _ 2 ' 



Thus for q < lo/c and q > lu/c the correlation functions are determined by the first and the second terms in Eq. ()19p . 
respectively. 



III. CALCULATION OF THE FRICTION FORCE 



The force which acts on the surface of body 1 can be calculated from the Maxwell stress tensor aij , evaluated at 
z = 0: 



1 

47r7o 



duj 



<Pq 



< E,E* > + < E*Ej > + < BiB* > + < B*Bj > 



-%(<E-E* > + <B-B* >) 



z=0 



Using Eqs. pi4p for the x - component of the force we get 

1 poo r j2 

11,1 d q q. 



Att 



duj 



(27r)2 fc2 



[{k..+kl){{\w,\') + {\w.\') 



(20) 



-{\vp P) - (I vs P)) + [k, - kl) {wpv; + w,vt - c.c)] 



(21) 



Substituting Eqs. (jlSlllSp for the amplitudes of the electromagnetic field in Eq. (f2T|) . and performing averaging over 
the fiuctuating electromagnetic field with the help of Eq. P^ . we get the a:-component of the force 

F. = a,, = A r du; f d'q^ [{q^ - (3kq,)' + /S^klq^] 

O"" Jo Jq<uj/c 1^1 



x[(g2 _ pkq^fil- I Rip n(l- I R' \^)\D,,\ 



+PXql{l- I Rip n(l- I i?L \')\Dsp\' + {P^ s)] {n^{J) - mic)) 



+ 2^^ r '^" />. c ^' W'^'' ~ ^^'^^' ^ /3^fc.^9^]e-2|'=^l'^ 




log V (m/s) 



FIG. 2: The velocity dependence of the frictional stress acting between two semi-infinite bodies at d = lOnm (1) and d = lOOnm 



(2), with parameters chosen to correspond to copper (r 
10) 



2.5 • 10" s 



1.6- 10'" s" 



(The base of the logarithm is 



: [{q^ - pkq^flmRiplmR'^plDssl^ ~ ^^klqllmRiplmR'^jD 



sp\ 



(pOs)](n2(w')-"i(^)). 



(22) 



where 



n.i{uj) 



1 



ahLo/kBTi _ I ' 



where Ti and T2 are the temperatures for bodies 1 and 2, respectively. The symbol {p O s) denotes the terms which 
can be obtained from the preceding terms by permutation of the indexes p and s. The first term in Eq. (j22p represents 
the contribution to the friction from propagating waves {q < w/c), and the second term from the evanescent waves 
{q > oj/c). If in Eq. ((22|) one neglects the terms of the order /3^ then the contributions from waves with p- and 
s- polarization will be separated. In this case Eq. ([22) is reduced to the formula obtained in^. Thus, to the order 
/3^ the mixing of waves with different polarization can be neglected, what agrees with the results obtained inS. At 
T = K the propagating waves do not contribute to friction but the contribution from evanescent waves is not equal 
to zero. Taking into account that n(— w) = — 1 — n{u)) from Eq. (|22p we get friction mediated by the evanescent 
electromagnetic waves at zero temperature (in literature this type of friction is named as quantum friction^) 



dqy 



dqx 



did 



qx 



[(?' - Pkqx) 



-P^klql]e 



■2\k^\ 



X [{q^ ~ l3kqxYlmRiplinR'2p\DssV " /S'fcf g;imi?iplmi?^ji?. 



{p O s)] 



(23) 



Fig. 2 shows the dependence of the frictional stress between semi-infinite bodies on the velocity V at different sepa- 
ration d. In the calculations the Fresnel formulas for the reflection amplitude were used with the Drude permittivity 
e for copper. The frictional stress initially increases with velocity, reaches a maximum at T^ ~ dksT /Ti^ and then 
decreases at large values of the velocity. 

The friction force acting on a small particle moving in parallel to a flat surface can be obtained from the friction 
between two semi-infinite bodies in the limit when one of the bodies is sufficiently rarefied. For d <C At = cfi/ksT, 
in Eq. (j22p we can neglect by the first term and in the second term we can integrate over the whole q-plane, and 
put kz ~ iq. We will assume that the rarefied body consists of small metal particles which have the electric dipole 
moment and the magnetic moment. The dielectric permittivity and magnetic permeability of this body, say body 2, 
is close to the unity, i.e. £2 ~ 1 ~^ AirnaE <?C 1 and /i2 — 1 — ^ AirnaH <?C 1, where n is the concentration of particles 
in body 2, as and an are their electric and magnetic susceptibilities. To linear order in the concentration n the 
reflection amplitudes are 



R2p — 



£2kz 



y/e2P-2k'^ - q^ ^ £2 - 1 



£2kz 



y/e2p.2k'^ - cp- £2 + 1 



lixna^ 




log V (m/s) 

FIG. 3: The velocity dependence of the friction force acting between smaU copper particle with radius R = lOnm moving above 
a copper sample at the separation d = 20nm. The contributions from the electric dipole moment and the magnetic moment 
are shown separately. (The base of the logarithm is 10) 



i?2 



fi2kz - \j£2li2k'^ - q^ 



M2 



^^2kz + Ve2A^2fc^^V ^^2 + 1 



2'KnaH 



The friction force acting on a particle moving in parallel to a plane surface can be obtained as the ratio between the 
change of the frictional shear stress between two surfaces after displacement of body 2 by small distance dz, and the 
number of the particles in a slab with thickness dz: 



ppart 



da\\{z) 



ndz 



= ^r / duj I d^q- 



qxq 



PW, 



e-29rf(^2(^')_„^(^)) 



\q (Imi?plma^ + Imi^^Ima'jLf) + /3 qyilvciRpluvau + hnRshi\a^)\ 



(24) 



where a', 



E(H) 



bj^ 



'^E(H){'^')- For a- spherical particle with radius R the electric and magnetic susceptibilities are given 



(25) 



aH 





t_fc, 


^ £ + 2 


f?3 

2 


1- 


-i?2.2+i?.^°*H 



(26) 



where k = k^/e. Fig. 3 shows the velocity-dependence of the friction force which acts on a small copper particle 
with R — lOnm, moving above copper sample at d = 20nm. The contributions from the electric dipole and magnetic 
moments are shown separately. At small velocities the contribution from the magnetic moment is seven orders of 
magnitude larger then the contribution from the electric dipole moment. 

IV. CALCULATION OF THE CONSERVATIVE FORCE BETWEEN MOVING BODIES 



The z-componcnt of the force is determined by zz-componcnt of the Maxwell stress tensor: 

I r . [■ d^q k. 



A-K 



duj 



(27r)2 A:^ 



[{kz + k:){{\w,\^) + {\wsn 



+ (I Vp n + (I Vs P)) + (fc, - kl) {{wXp + WsV; + C.C)] 



(27) 



One has to subtract the infinite vacuum contribution to the force which does not depend on separation d^^. Substi- 
tuting (jlSmSp into Eq. (|27|) , and averaging over the fluctuating electromagnetic field with the help of Eq. (fT9| , and 
after subtraction of the vacuum term, we get the z-component of the force: 

fi r°^ r k r 

+RisR'2sDpp] - ^^klqllRipR'^^Dsp + RisR2pDps]}[l + ni(cj) + n2(cj')] 



x{{q' - /3fcg,)2[(l- I R,p n(l+ I i?^p \^)\D,s\' - (1 O 2,c. o io')] 
+f3^kX[{l- I Rip |2)(1+ I R',^ \^)\Dsp\^ - (1 o 2,c. f^ u:')] + {p ^ s)} {n,{u) ~ n^^J)) 

x{(g2 _ ^fcg^)2[Inii?ipRei?^p|D,,|2 - (1 o 2,tj f^ J)\ ~ f? klq^y\iYiiRxp^QR!^^\D sp"^ - (1 O 2,tj O a;')] 

+ (pOs)}(ni(Lj)-n2(w'))- (28) 

At Ti = T2 = K, Eq. dMl) takes the form 

{(g2 - pkq,f[R^pR'^^D,, + R^.R'^^Dpp] 
-P k^qy[RipR'2sDsp + RisR'2pDps]> 

^ ^°° dqy r dq, r'' dJ^[{q^ ~ Pkq^f + P^klql]e~^\^'\^ 



47r3j_^ JO JO 

x{(g2 _ pkq,f[lxnRipReR'^^\Dss? - (1 O 2,lj o w')] - /3'fc^g^[Imi?ipRei?iJD,p|2 - (1 o 2,cj o u')] 

+ {p^s)}{ni{uj)~n2{uj')). (29) 

If in Eq. ([25)1 one neglects the terms of the order /3^ , then, as in the case of friction, the contributions from the waves 
with p- and s- polarization will be separated. In this case Eq. (j28p reduces to 

h f°° f f 1 

F^ = ^Re / dw / d^qk. 



47r3 Jo J ^ '\R^^'R',-'e-^^^^d-l 

[1 + ni{uj) + n2{uj')] 



Ru'KS'e-^^'^'^ - 1 



n r d^ ( ^2„,. r [(1-1 ^ipr)(i+ 1 ^2pr)- (1^2,0.^^0.')] 



^9fc.{- 



167r3 7o Jq<u:/c ^ \Dpp\2 

+ (pO s)|(7ii(a;)-?i2(w')) 



47^ Jo J9>a;/c 



+ (pOs) WniH-n2(^')) (30) 



[Imi?ipRei?^p ~ (1 o 2, w o w')] 

If we put V = Q and use the Fresnel's formulas for the reflection amplitudes, then Eq. (l30t reduces to the formula 
obtained by Lifshita^. Lifshitz have shown that at Ti = r2 = K it is convenient to transform w-integration along 
the real axis into the integral along the imaginary axis in the upper half of the complex cj-plane. For a rarefied body, 
similarly as in Sec. IIIIl for d ^ Ti/ksT from (pS)) we get the van der Waals interaction between a small particle and 
plane surface: 



h f°° f q^e '^ r 

Fz = ^Im / duj I d^q— — ^-^ q^lRpOiE 

71"^ Jo J q^ - P^qti •- 

+i?,a^] + P^qllRpa'fj + i?sa^]}[l + n^{uj) + 712(0;')] 



'2^ q ^-2qd 





.2 ,_ du^jd q^^^e 



x{q [ImRpKea'^ — ReJ^pIma'^] + (3 q^, [Imi?pRea^ — Imi?sRca^] 

+ (p ^ s, E ^ H)} imiuj) - n2{Lj')) . (31) 

V. CALCULATION OF THE RADIATIVE HEAT TRANSFER BETWEEN MOVING BODIES 

The radiative energy transfer between the bodies is determined by the ensemble average of the Poynting's vector. 
In the case of two plane parallel surfaces the heat flux across the surface 1 is given byi^ 

(Si.(r))^ = (c/8^)(E(r)xB*(r))„ + c.c. 

= ^((E(r).-fE*(r))-c.c| . (32) 

Using Eqs. pi4p and ([5^ we get equations for the heat obtained by the body 1 which are very similar to Eq. (PT|) 
and ^■. 

- (I vp P) - (I V, P)) + {K - K) {{w,v; + w^vl - c.c)\ . (33) 

After averaging of the product of the components of the fluctuating electromagnetic field in the same way as in Sec. 
dnil we get 



x[{q^-Pkq^f{l~\R^,\^){l-\R'^^\^)\D,,\'' 
+P^klql{l- I i?ip n(l- I R'^^ \^)\Dsp\^ + {p^ s)] (n2(w') - ni(L.)) 

x[{q^ - pkq^flTnRiplTnR'^^\Dss\'' - l3^klqllmRipImR'^,\Dsp\^ 

+ {p^s)]{n2{uj')~ni{ij)). (34) 

Eq. (|M)) generalizes the equations for the heat transfer between two surfaces which are at the rest in the X-reference 
franie^'^^, to the case when the surfaces are moving relative to each other. There is also the heat 5*2 obtained by the 
body 2 in the X'-reference frame. Actually, Si and 5*2 are the same quantities, looked at from different coordinate 
systems. These quantities are related by the equation: 

F^V = 5i + S2/1. (35) 

For the limiting case of rarefied body, similarly as in Sec. IIIIl from ([34]) we get the contribution from evanescent 
waves to the heat adsorbed by a semi-infinite solid at d ^ At in the iiT-reference frame: 

^=^ Tdc. fd?q^^^e'^^^[n,{u')-ni{u:)) 

TT^ Jo J r - P % 

X [g2(lmi?plma^ + Imi^Jma'j:^) + /3^(j^(lmi?plma'^ + Imi?Jma^)] (36) 

The heat adsorbed by a particle can be determined from Eq. psp . 



VI. CALCULATION OF THE FRICTION FORCE ON A SMALL NEUTRAL PARTICLE MOVING 

RELATIVE TO BLACK BODY RADIATION 

We consider a small neutral particle moving relative to black body radiation. We introduce two reference frame K 
and K' . The thermal radiation is in equilibrium in the X- reference frame and the particle is at rest in the if'-reference 
frames. We assume that the particle moves with velocity V along the z-axis. The relation between the cc-components 
of the momentum in the different reference frames is given by 

p^ = {p', + f3Eo/c)^ (37) 

where Eq is the rest energy of the particle. The rest energy can change due to thermal radiation of the particle. From 
Eq. dSZl) we get 

dpx _dp'V dEo . . 

dt ~ dt' ^ c2 dt' ^^^^ 

According to the Einstein law 

dmo _ 1 dEo 
dt' ~ c2 dt' ^'^^' 



where mo is the rest mast of the particle. Taking into account that 



from Eqs. pSMOp we get 



dp^ d{moV-f) dnio d{Vj) 

— ;— = ; = , , V + mo ; (40) 

dt dt dt' dt ^ ' 



d{Vl) dp'. 



In the rest reference frame, due to symmetry, the total radiated momentum from the dipole and magneto-dipole 
radiation is identically zero. Thus, the change of momentum of the particle in the rest reference frame is determined 
by the Lorenz force F', acting on the particle from the external electromagnetic field associated with the thermal 
radiation observed in this reference frame. The dynamics of the particle in the iiT-reference frame is determined by 
the equation 

-0^ = ^^ (42) 



Eq. (l42t does not contain force from the thermal electromagnetic field radiated by the particle. Thus, from Eq. p2)) 
it follows that, contrary to the claim of the authors of^'', the thermal radiation of the particle can not produce any 
acceleration. In the iiT'-reference frame the Lorenz force on the particle is determined by the expressio n^^i^^ 

K = ^(Pe • E'*(r'))r'=ri + ^ (Pm ' B'* (r')>r'=ri (43) 

Writing the electromagnetic field as a Fourier integral, and taking into account that 

p^=aB(c.')E'(r;,)e— '*', 



P™ = aff (^')B'(r;,)e 



-iuj' t' 



we get 






L^J Wf^'-" '"""^"'^ ^ """^"'^^ ^^' ■ ^'*^"''' 



(44) 



where we have taken into account that for plane waves (E' • E'*)^'k' = (B' • B'*)tj'k'- When we change from the K'- 
reference frame to the iC-reference frame (E' • E'*)i^/k' is transformed as the energy density of a plane electromagnetic 
field. From the law of transformation of the energy density of a plane electromagnetic field?- we get 

(E' • E'*)„,k' = (E • E*)„k f^") . (45) 

According to the theory of the fluctuating electromagnetic field'^^ 

(E • E*)„k - 47r2nfc {6{- - fc) - <5(- + fc)| [1 + 2n(a;)] (46) 

L c c J 

Taking into account the invariance of the square of the four-wave vector (w/c)^ — fc^ = (w'/c)^ — k'^ Eq. (|151) can be 
rewritten in the form 

(E . E*)^k = ^^^ |^(^ - k') - 5(^ + k')\ [1 + 2n(c^)] (47) 

Substitution Eqs. (|45l47p in Eq. pi)) and integration over a;' gives 

he f 
Fx = -p-—^ I d^kkkx [Ima£;(cfc) +linaH{ck)] 



X {n[7(cfc + Vk,)] - n[-f{ck - Vk,)]} , (48) 

where we have omitted index of prime and have taken into account that uj = {uj' + k'^V)j. Introducing the new 
variable lo = ck, (|48p can be transformed to the form 

2h r°° , , /■"/ 



Fx == — J / duju'^ / dk^^kx [Ima£;(a;) + Imaij (w)] 
""c Jo Jo 



X {n[7(cfc + Vkx)] - n[7(cfc - Vk,)]} (49) 

Eq. P^)) generahzes the resuh obtained in^^ to the case of large velocities, and includes the contribution from 
the magnetic moment. For metallic particles the contribution from the magnetic moment exceeds substantially the 
contribution from the electric dipole moment. At small velocities F^ ~ ~rV^, where 



r = ^JZ3 / '^\^]'^'^ P™"^ ^'^^ + ^"""^ '^'^)] ^^^"^ 



Sttc^ 



For metals with Attu ^ ksT /% and for c^/2TiakBT ^ R, where a is the conductivity, from Eqs. psp and (|26|) we get 

Ima£;(w) w i?^-^ (51) 

47rCT 

Imaniuj) = ^^^^ (52) 

Setting the friction coefficient T to mo/r, where r the relaxation time, and using ttiq = 47ri?^p/3, from Eqs. ([501- 152)1 
we get 

, , hR <jR , , 

where r^^ and r~^ are the contributions to the friction from the electric dipole and magnetic moments, respectively. 
For T = 300 K, p w lO"* kg/m^, a sa 10^^ s'^ from Eqs. (|51I54I) we get t^ - lO^^ s and -r,„ - lO^^ s. When the 
conductivity decreases Te also decreases and reaches minimum at 2'Ka w ksT/h. At T — 3000 K this minimum 
corresponds to about a day (r™™ « 10^ s). InSS the same relaxation time was obtained for Ba^. 

VII. CONCLUSION 

In this paper within the framework of unified approach we have calculated the Casimir-Lifshitz interaction, the van 
der Waals friction force and the radiative heat transfer at nonequilibrium conditions, when the interacting bodies are 
at different temperatures, and they move relative to each other with the arbitrary velocity V. In comparison with the 
existing literature we have studied more general nonequilibrium conditions. Our study was focused on the surface- 
surface and surface-particle configuration. We have found the exact solution of problem about the determination of 
the fluctuating electromagnetic field in the vacuum gap between two fiat parallel surfaces moving relative to each 
other. Knowing the electromagnetic field we have calculated the Maxwell stress tensor and the Poynting vector which 
determine the friction and conservative forces, and the heat transfer between the solids, respectively. For the heat 
transfer and the conservative force our treatment generalizes the results obtained for bodies at rest to the case of 
bodies which move relative to each other. We have found that the velocity dependence of the considered phenomena 
is determined by Doppler shift and can be strong for resonant photon tunneling between the surface modes. This 
effect can be used for the precision determination of energy of the surface modes. We have shown that relativistic 
effects produce a mixing of the s and p-wave contributions to the forces and the heat transfer. This relativistic effect 
is of the order {V/c)"^. If one neglects by terms of order {V/c)"^, the different polarizations will contribute to the 
forces and the heat transfer separately. The velocity dependence of the conservative force is much weaker than the 
velocity dependence of the friction force and the heat transfer. For the conservative force the important range of the 



frequencies is determined by the plasma frequency ujp which is much larger than the Doppler shift Aw = QxV ~ V/d 
for practically all separations and velocities. Therefore the velocity dependent component of the Casimir-Lifshitz 
force will be considerably smaller than the V = component. The same is true for the thermal component of the 
Casimir-Lifshitz force. The thermal component was measured recently in—. Thus, one can expect that the velocity 
dependent component also can be measured. For the friction force and the heat transfer the important range of the 
frequencies is determined by the thermal frequency wt = ksT/h which can be of the same order as the Doppler shift. 
The friction force initially increases when the velocity increases, reaches maximum at ^ ~ wt, and then decreases. 
From the limit when one of the bodies is rarefied, we have calculated the interaction force and the heat transfer 
for a small particle outside a flat surface. For a particle we have taken into account the contribution to the forces 
and the heat transfer from the electric dipole and magnetic moments. For metallic particles we have found that the 
contribution to the friction from the magnetic moment exceeds the contribution from the electric dipole moment by 
several orders of magnitude. We have presented a relativistic theory of the friction force acting on a particle moving 
relative to the black body radiation. We have shown that the thermal radiation of the particle does not produce any 
acceleration and that the friction force acting on the particle is determined by external electromagnetic field associated 
with the black body radiation. For a metallic particle moving relative to the black body radiation the friction force 
increases initially when the conductivity decreases, reaches maximum at 27r cr « cut, and then decreases. 
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